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Confinement may be more easily demonstrated at finite temperature using the 
Polyakov loop than at zero temperature using the Wilson loop. A natural mechanism 
for confinement can arise via the coupling of the adjoint Polyakov loop to F^^,. We 
demonstrate this mechanism with a one-loop calculation of the effective potential 
for SU (2) gluons in a background field consisting of a non-zero color magnetic field 
and a non-trivial Polyakov loop. The color magnetic field drives the Polyakov loop 
to non-trivial behavior, and the Polyakov loop can remove the well-known tachyonic 
mode associated with the Saviddy vacuum. Minimizing the real part of the effective 
potential leads to confinement, as determined by the Polyakov loop. Unfortunately, 
we cannot arrange for simultaneous stability and confinement for this simple class of 
field configurations. We show for a large class of abelian background fields that at 
one loop tachyonic modes are necessary for confinement. 
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I. INTRODUCTION 



While the problem of quark confinement has been with us for some time [0 , a completely 
satisfactory explanation of the mechanism has been elusive. It is widely held that unbroken 
non-Abelian gauge theories functions as dual superconductors, with some analog of the 
Meissner effect in type II superconductors confining chromoelectric fiux into fiux tubes. 
A fully satisfactory theory of quark confinement would be Lorentz invariant and gauge 
invariant, and have correct renormalization group behavior, which suggests that a successful 
theory must maintain a tight connection with perturbation theory. 

Confinement is easiest to demonstrate at finite temperature, where it can be studied using 
the trace of the Polyakov loop V = T exp i (HAq . In pure gauge theories, without fun- 
damental representation quarks, confinement is equivalent to the vanishing of the expected 
value of the trace of the Polyakov loop in the fundamental representation, TrpV = 0. In 
contrast, demonstrating confinement at zero temperature requires showing that the Wilson 
loop has area law behavior for large areas. To use a simple analogy, it is much easier to 
study one-point functions than two-point functions. In pure SU{N) gauge theories, there is 
a global Z{N) symmetry, which, if unbroken, enforces the vanishing of the Polyakov loop. 
The spontaneous breaking of this Z{N) symmetry at high temperature is observed in lattice 
gauge theory simulations. In accord with these results, a successful theory of confinement 
should also show that the expected value of the trace of the Polyakov loop vanishes at low 
temperatures, and is non-zero at high temperatures. 

The perturbative effective potential for the Polyakov loop, or equivalently, Aq in an 
appropriate gauge does not indicate confinement @J^. The general result for SU{N) is a 
quartic polynomial in Aq. The minimum of this effective potential is always at Aq = 0, 
corresponding to trpiV) = N; in the pure gauge theory without quarks, there are also 
equivalent solutions related by the global Z{N) symmetry. There is a simple physical reason 
behind this: when Aq = 0, the free energy associated with the quark-gluon plasma is 
minimized. 
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One possible origin of confinement is that at low temperatures different regions of space 
are in different phases, in such a way that trp{V) averages to zero. This possibility 
is realized in the strong-coupling region of Z{N) lattice gauge theories, where individual 
Polyakov loops can only take on values in Z{N), but the average over an infinite volume is 
zero. This Z{N) confining phase is connected to a phase of SU{N) lattice gauge theories: if 
the Wilson action PfSf has added to it an adjoint term PaSa, then the Z{N) gauge theory 
can be obtained in the limit (3a ^ oo with f^p fixed. The region Pa large and (3f small 
realizes confinement in the same manner as the Z{N) model. However, this region of the 
Pf-Pa plane is separated from the continuum limit fixed point by a line of first-order phase 
transitions [Q. In contrast, examination of lattice field configurations generated by Monte 
Carlo simulation along the Pa = ^ axis indicate that the magnitude of individual Polyakov 
loops is small at low temperatures and increases dramatically at the confinement transition. 
This argues against the possibility that confinement occurs because different regions of space 
are in different Z]y phases. 

We advocate another possibility: that the coupling between the local gauge field F^^^, and 
the adjoint Polyakov loop produces an effective action which leads to two different phases. In 
the low temperature phase, there is a magnetic condensate and the Polyakov loop indicates 
confinement; in the high temperature phase the magnetic condensate vanishes, the Polyakov 
loop indicates deconfinement, and the contribution of the thermal gauge boson gas dominates 
the free energy. This occurs because finite temperature effects from gluons naturally couple 
the local gauge field to adjoint representation Polyakov loops 

In order to produce confinement, the low temperature phase must minimize its free 
energy by minimizing the expected value of the trace of the adjoint Polyakov loop. For 
SU{N), the fundamental and adjoint representation Polyakov loops are related by 



Clearly, when TtaIV) assumes its minimum value of —1, Trp{V) assumes the value 0. As 
seen in previous calculations, at sufficiently high temperature, contributions from thermal 



TrA{V) = \TrF{V)\^-l. 
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gluons will dominate the free energy, driving the Polyakov loop towards its maximum value. 
In our scenario, this also drives the magnetic condensate to zero. The coupling of the 
gauge field to the adjoint Polyakov loop naturally incorporates the behavior known from 
lattice simulations: both the Polyakov loop and plaquette expectation values change at the 
deconfining phase transition of SU{N) lattice gauge theories. 

The idea that more than one field must be considered to describe the behavior of finite 
temperature QCD has previously been proposed by Campbell, Ellis and Olive [Il5| , [l6| . In 



their effective theory, a phenomenological coupling of the gluon condensate to the chiral 
condensate was postulated. However, this approch does not include the important role 
played by the Polyalov loop in finite temperature chiral behavior. It is easy to show that 
there is a coupling between the chiral condensate and Polyakov loops. For example in 
references and , we showed using lattice perturbation theory the role that the Polyakov 
loop plays in finite temperature corrections to the renormalization of the gauge field coupling 
constant. Differentiation of our result with respect to the quark mass shows perturbatively 
that the finite temperature corrections to the chiral condensate depend on the Polyakov loop 
as well as the gluon condensate. We have also shown in detail the interplay between the 
Polyakov loop and the chiral condensate for a variant of the Nambu-Jona-Lasinio model. 
0J^. A truly satisfactory treatment of the critical behavior of finite temperature QCD, 
including the effects of quarks, is likely to require treatment of the gluon condensate, the 
chiral condensate and the Polyakov loop. 

II. CALCULATION OF THE FREE ENERGY 

The observation of Saviddy that the perturbative QCD vacuum is unstable with 
respect to the formation of a constant chromomagnetic field is one proposed origin for 



a magnetic condensate in QCD. However, Nielsen and Olesen |T0[ showed that such field 
configurations are themselves unstable, due to a tachyonic mode in the one- loop determinant. 
Several authors |]TT]-[T^ have extended the work of Nielsen and Olesen to finite temperature, 
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treating the tachonic modes in various ways, including simply ignoring them. 

We demonstrate the main points of our proposed mechanism by incorporating the effects 
of a non-trivival Polyakov loop, treating all modes exactly at one loop. The tachyonic mode 
can be removed by a non-trivial Polyakov loop. The free energy of an SU(2) gluon gas moving 
in a background field consisting of a constant color magnetic field and and a constant non- 
trivial Polyakov loop is obtained by calculating the one-loop, finite temperature effective 
potential in the background field gauge; see fl^ and for pedagogical introductions to 



the background field method. This is equivalent to summing up three contributions to the 
free energy: the classical action, the zero-point energies, and the free energy of the gluon 
gas. All finite temperature effects reside in the last term. As we have shown elsewhere, 
Polyakov loops appear naturally at one loop in an image expansion of finite temperature 
determinants ||5|-Q. Here, similar results are achieved by directly computing the functional 
determinant. 

We choose the color magnetic field and the Polyakov loop to be simultaneously diagonal, 
and we let the color magnetic field H point in the X3 direction. We take the Polyakov loop 
to be specified by a constant Aq field, given in the adjoint representation by 

A„ = |r3. (2.1) 

The trace of the Polyakov loop is then given by 

Trj.(P) = 2 cos(0/2) (2.2) 
in the fundamental representation and by 

TrA(P) = 1 + 2cos((/>) (2.3) 
in the adjoint representation. The external magnetic field we take to have the form 

A2 = ^HxiT3 (2.4) 
which gives rise to a chromomagnetic field 
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2 ^ 



(2.5) 



As explained, for example, in [ITU |, the external field gives rise to Landau levels in the gluon 
functional determinant. 

The one-loop contribution to the free energy has the form [|lO| , |ri[] : 



2xi£ J:^^ f*ln 



^n-|j +2H(m + ^±l]+k:, 



+ E 



1 



In (ul + P 



(2.6) 



(3 J (27r)3 

where the a;„ = 2Tm/ (3 are the usual Matsubara frequencies, and where the terms 2H{m + 
1/2 ± 1) are the allowed Landau levels of the gauge field with its spin coupled to H. The 
explicit factor of 2 in front of the first term in the sum results from the trace over color 
degrees of freedom. 

When = 0, the n = and m = modes give rise to tachyonic modes for k^ sufficiently 
small; these in turn give rise to an imaginary part in the free energy. We observe that these 
same modes will give a strictly real factor to the determinant provided 

> P^/H. (2.7) 

However, if is too large, it may cause the n = 1 mode to become unstable. Hence, the 
general criterion for stability is 



(3Vh < < 277 - pVh 



(2.8) 



A standard product representation 



n 

k 



1 + 



V2A;7r 



cosh{x) — cos (a) 



1 — cos(a) 



allows us to write the effective potential in the form: 

H r dk 



1 1 °° 



d^k 



, 27rP J 2vr 



In {cosh [/3ui±{m^ k)] — cos(/90)} 



+ 



(3 J (27r)3 



In {cosh [p\k\\ - l} 



(2.9) 



(2.10) 
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where 



ulim, k) = 2H (m + ^ ± l) + 



(2.11) 



which can be negative. Note that the classical free energy of the magnetic field, H'^/2g'^, is 
included in this expression. Expanding the logarithm and discarding an irrelevant constant, 
the contribution of the Landau levels to the the effective potential can also be written in 
the form: 



-EE- 

2 97T 
^ m=0 c 



dk 
2^ 



2 ^ cosfnf/)^) 



-nfiui± (m,k) 



n=l 



n 



(2.12) 



After performing the sum over m in the equation above, we regulate the divergent por- 



tions of the resulting integral using the identity |T0 



r(/i - v) 

This yields a renormalized effective potential with real component 



(2.13) 



IIH^ (H 



{Hfl^ g cos(n: 



and imaginary component 



cosini 

n=l 



n 



E 



Ki{nf3VH)--Yi{nl3VH) 

27r2 



n 



V2m + 3Ki[nf3^{2m + 3)H] 

m=0 



90/3^ 



(2.14) 



tt2 f7'3/2 oo 

Vj = -— -^y 



co^[nc 



n 



Ji{nf3VH) 



(2.15) 



where is defined in Eq. (|2.1| ). 

It is far from obvious that the potential we have found is real when the constraint of 
Eq. ( |2.8| ) is satified. We have checked numerically that the imaginary part of the potential 
behaves as it should. In figure 1, we show the imaginary part for \/H / ^ = 0.5 and T//i = 
0.25. The theory is stable in the region given by Eq. ( |2.8| ). 

A similar calculation has been performed for different reasons by Starinets, Vshivtsev, 
and Zhukovskii [^. Our results are in disagreement, with Bessel functions interchanged 
between the real and imaginary part. However, our results are in numerical agreement with 



the exact result for Vj derived by Cabo, Kalashnikov, and Shabad p2[ for the case = 0, 
and our result for Vj is zero for the range of (p determined by Eq. (|2.8| ). We are confident 
our results are correct. 

Note that we have carefully chosen the overall additive constant of the free energy such 
that V = when = and T = 0. 

III. MINIMIZATION OF RE(V) 

We are interested in showing that the system can lower its free energy by having = vr, 
or equivalently TrpP = 0. The simplest scheme is to minimize Re{V) over all allowed values 
of (f) and H, ignoring the imaginary part of V. If the global minimum of Re{V) gives rise to 
an imaginary part, that will indicate that this field configuration is unstable, however. 

Some insight into this issue can be obtained by noting that for very low temperatures, the 
finite temperature contributions to Re{V) will be dominated by the Yi{nl3\/H) terms. Using 
the asymptotic form for large argument, we see that they will favor the adjoint Polyakov 
loop being at its minimum, and hence favor confinement. 

We have examined this issue numerically and found that at low temperatures, minimiza- 
tion of Re{y) leads to = vr being preferred. This is shown in figure 2, which plots the 
real part of the effective potential versus -ff/yU^ at T//i = 0.25 for both = and = vr. 
Figure 3 plots V/-, the imaginary part of the effective potential for the same parameters. 
Unfortunately, examination of V/ shows that the lowest minima is not stable, lying just to 
the right of the stable region. 

The Saviddy vacuum is preferred over the perturbative vacuum only for sufficiently low 
temperatures. Figure 4 shows the behavior of V at T//i = 0.73 and = vr compared with 
Vqti the one-loop effective potential for a gluon gas at = and = 0. Above the 
pseudocritical temperature at approximately Tc/fi = 0.73, the perturbative solution with 
H = and = has lower free energy. Naively, this would indicate a first-order phase 
transition. However, the magnitude of H/ fj,"^ is so large that a one-loop calculation cannot be 



8 



considered reliable. Further, examination of Vj shows that the constant field configuration 
which minimizes V is always unstable. All that can be said is that at low temperatures, the 
perturbative solution is definitely not the state of lowest free energy. 

It is amusing to note that the behavior at very low temperatures reveals the existence 
of a familiar phenomenon associated with Landau levels, the de Haas- van Alphen effect. 
Fig. 5 plots the free energy as a function of H/fi^ at T/fi = 0.1 for both = and 
(f) = IT. The minima are the signs of the characteristic oscillatory components of the magnetic 
susceptibility. 

Because of the possibility of fixing the external field H in lattice simulations, it may be 
of interest to examine the analytical behavior at fixed H. At fixed H, = vr is not necessarily 
the minimum value of V. This is demonstrated clearly in Fig. 6 and Fig. 7, which show the 
V and Vi as a function of at T/fi = 0.25 for H/j^ = 0.40,0.65 and 0.90. For H/jj,'^ = 0.40, 
there is a region of where Vi vanishes. The minimum of V lies at or very near the boundary 
where Vj develops an imaginary part. This particular case appears similar to the behavior 
found by van Baal, who examined self-dual solutions on a hypercube in the presence of a 
non-trivial Polyakov loop p5|] . 



As we have seen above, the terms in the free energy which favor confinement come from 
the contribution of the would-be tachyonic modes to the free energy. We can generalize this 
result slightly. Consider an arbitrary time-independent vector potential of the form 



IV. GENERALIZATIONS 



1(f) rs 



(4.1) 



We write the eigenvalues that enter into the functional determinant in the form 




(4.2) 



where lo"^ can be positive or negative. 
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We write the positive eigenvalues as e^, and the negative eigenvalues as — x^- Following 
the identical reasoning used above, the positive eigenvalues contribute to the free energy a 
term 



(4.3) 



where p+ (e) is the eigenvalue density. It it easy to see that such eigenvalues always favor 
= 0. On the other hand, the negative eigenvalues contribute to the free energy a term 

2 



X 



= C + Jdx p{x) In [cos^ {(3x) - cos^ (0)] , (4.4) 

and such terms may favor 7^ 0. Thus in this simple generalization we have an indication 
of the likely importance of the negative modes in bringing about confinement. 

It is interesting to note that the hard thermal loop resummation of Braaten and Pisarski 
2^ drops precisely those features which we are advocating as the origin of confinement at 



low temperatures. The tachyonic modes in the Saviddy vacuum arise at the one-loop level 
from the spin coupling of the gauge field to the external magnetic field. However, the term 
accounting for this spin coupling is discarded in the hard thermal loop approximation, since 
the term is manifestly soft. 



V. CONCLUSIONS 

We have seen that finite temperature gauge theory is the best place to prove confinement, 
because one need only show that the expected value of the Polyakov loop vanishes. Models 
and mechanisms for which this cannot be demonstrated are not good candidate explanations 
of confinement. We have demonstrated these considerations with a one- loop calculation of 
the effective potential for SU(2) gluons in a special background field consisting of a non- 
zero color magnetic field and a non-trivial Polyakov loop. The color magnetic field drives 
the Polyakov loop to non-trivial behavior, and the Polyakov loop can remove the well- 
known tachyonic mode associated with the Saviddy vacuum. Minimizing the real part of the 
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effective potential leads to confinement, as determined by the Polyakov loop. Unfortunately, 
we cannot arrange for simultaneous stability and confinement in this simple model. We are 
able to show for a class of abelian field configurations that tachyonic modes are necessary 
for confinement. 

The mechanism we have described here for pure gauge theories is also applicable when 
the effects of quarks are included. The quark determinant has several effects, the most 
important of which is the explicit breaking of Z{N) symmetry. Pure SU{N) gauge theories 
have an exact global Z{N) symmetry which, if unbroken, requires that the expectation value 
of TrpP vanish. Fundamental representation particles, such as quarks, explicitly break this 
symmetry. However, explicit perturbative calculations shows that the symmetry breaking 
terms are numerically small, even when the current quark mass is zero. Strong-coupling 
arguments ||26| suggest that in the low temperature phase, it is the constituent quark mass 
rather than the current mass which controls the strength of the Z{N) symmetry breaking, 
which would further reduce the effect. 

Other effects of quarks enter from the change in renormalization group beta function 
coefficients, and from the increase in the number of degrees of freedom in the high tempera- 
ture phase. At low temperatures, we expect that the competition between the gluonic terms 
which drive the trace of the Polyakov loop towards zero and the quark terms which drive it 
away from zero will produce a low temperature region where the Polyakov loop is small, in 
accord with lattice simulations. 

It is obvious that the ground state of QCD is not well modeled by a constant chromomag- 
netic field. However, we can extend the insight we have gained into the possible origins of 
confinement by examining the confining properties of other background field configurations. 
There are several interesting classes of configurations, including constant field potentials. 
If the different space time components do not commute, this leads to constant color fields 



which are inequivalent to the class studied here ||2^. Another interesting class of background 
field configurations is abelian monopole configurations. It would also be very interesting to 
know if a suitably chosen ensemble of random background field configurations confines. 

11 



ACKNO WLED CEMENTS 

We wish to thank the U.S. Department of Energy for financial support under grant 
number DE-FG02-91-ER40628. 



12 



REFERENCES 

[1] S. Mandelstam, Phys. Kept. 67, 109 (1980). 

[2] N. Weiss, Phys. Rev. D 24, 475 (1981). 

[3] N. Weiss, Phys. Rev. D 25, 2667 (1982). 

[4] G. Bhanot and M. Creutz, Phys. Rev. D 24, 3212 (1981). 

[5] P. N. Meisinger and M. C. Ogilvie, NucL Phys. B (Proc. SuppL) 42, 532 (1995). 

[6] P. N. Meisinger and M. C. Ogilvie, Phys. Rev. D 52, 3024 (1995). 

[7] P. N. Meisinger and M. C. Ogilvie, Nucl. Phys. B (Proc. Suppl.) 47, 519 (1996). 

[8] P. N. Meisinger and M. C. Ogilvie, Phys. Lett. B 379, 163 (1996). 

[9] G. K. Saviddy, Phys. Lett. 71B, lr363 (1977). 
[10] N. K. Nielsen and P. Olesen, Nucl. Phys. B144, 376 (1978). 
[11] M. Ninomiya and N. Sakai, Nucl. Phys. B190 [FS3], 316 (1981). 
[12] B. MuUer and J. Rafelski, Phys.Lett. lOlB, 111 (1981). 
[13] W. Dittrich and V. Schanbacher, Phys.Lett. lOOB, 415 (1981). 
[14] J. Kapusta, Nucl. Phys. B190 [FS3], 425 (1981). 

[15] B. A. Campbell, J. Ellis, and K. A. Olive, Nucl. Phys. B345, 57 (1990). 
[16] B. A. Campbell, J. Ellis, and K. A. Olive, Phys. Lett. B235, 325 (1990). 
[17] L. F. Abbott, Acta Phys. Pol. B13, 33 (1982). 

[18] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory, Addison- 
Wesley (1995). 

[19] W. Greiner and A. Schafer, Quantum Chromodynamics, Springer- Verlag, Berlin (1994). 

13 



[20] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Fifth Edi- 
tion, Academic Press, New York (1994). 

[21] A. O. Starinets, A. S. Vshivtsev, and V. Ch. Zhukovskii, Phys. Lett. B 322, 403 (1994). 

[22] A. Cabo, O. K. Kalashnikov, and A. E. Shabad, Nucl. Phys. B185, 473 (1981). 

[23] J. I. Kapusta, Nucl. Phys. B190 [FS3], 425 (1981). 

[24] E. Braaten and R. D. Pisarski, Nucl. Phys. B337, 569 (1990). 

[25] P. van Baal, Nucl. Phys. B (Proc. Suppl.) 47, 326 (1996). 

[26] A. Gocksch and M. Ogilvie, Phys. Rev. D 31, 877 (1985). 

[27] L. Brown and W. Weisberger, Nucl. Phys. B 157, 285 (1979). 



14 



FIGURES 

FIG. 1. Vi as a function of cf) for ijV2/^ = 0.5 and T/fi = 0. 
FIG. 2. V at T/iJ, = 0.25 for ^ = and phi = it. 
FIG. 3. Vi at T/n = 0.25 for cf) = and p/iz = tt. 
FIG. 4. y at T/n = 0.1 for (?!) = and p/tz = tt. 
FIG. 5. y at T//X = 0.73 for ^ = tt compared with Vqt- 
FIG. 6. F as a function of (f) at T//x = 0.25. 
FIG. 7. V/ as a function of cf) at T/n = 0.25. 
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Figure 1: Vj as a function of (|) for H^^^/ji = 0.5 and T/ji = 0.25 




Figure 2: V at T/|Li = 0.25 for (|) = 0,^ 




Figure 3: Vj at T/|Li = 0.25 for ^ = 0,7C 




Figure 4: V at T/|Li = 0.1 for (|) = 0,7C 




Figure 5: V at T^ill = 0.73 with ^ = 71 compared to V, 
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Figure 6: V as a function of (|) at T/[i = 0.25 
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Figure 7: Vj as a function of (|) at T/jj. = 0.25 



